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Reaction-diffusion equation for quark-hadron transition in heavy-ion collisions 
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Reaction-diffusion equations with suitable boundary conditions have special propagating solutions 
which very closely resemble the moving interfaces in a first order transition. We show that the 
dynamics of chiral order parameter for chiral symmetry breaking transition in heavy-ion collisions, 
with dissipative dynamics, is governed by one such equation, specifically, the Newell-Whitehead 
equation. Further, required boundary conditions are automatically satisfied due to the geometry of 
the collision. The chiral transition is, therefore, completed by a propagating interface, exactly as 
for a first order transition, even though the transition actually is a crossover for relativistic heavy- 
ion collisions. Same thing also happens when we consider the initial confinement-deconfinement 
transition with Polyakov loop order parameter. The resulting equation, again with dissipative 
dynamics, can then be identified with the reaction-diffusion equation known as the Fitzhugh-Nagumo 
equation which is used in population genetics. We discuss the implications of these results for heavy- 


ion collisions. We also discuss possible extensions 
PACS numbers: 68.35.Fx, 25.75.-q, 12.38.Mh, 64.60.-i 

I. INTRODUCTION 

Dynamics of quark hadron transition is one of the most 
important issues in relativistic heavy-ion collisions, as 
well as in the universe. Earlier it used to be believed 
that the quark-hadron transition is first order even at 
low chemical potential (as in the early universe). This 
gave a very important proposal by Witten [TJ about the 
possibility of formation of quark nuggets due to the con¬ 
centration of quarks by moving phase boundaries at the 
quark-hadron transition. Dynamics of first order transi¬ 
tion also had important implications for heavy-ion colli¬ 
sions [2]. Subsequently, lattice results showed that the 
quark-hadron transition is not first order, rather it is 
most likely a cross-over for low chemical potential. This 
cross-over is believed to govern the dynamics of transi¬ 
tion in relativistic heavy-ion collisions at high energies. 
(Though for lower energies, transition may become first 
order when the baryonic chemical potential is sufficiently 
large.) 

For the dynamics of the phase transition, the most im¬ 
portant difference between a first order transition and a 
cross-over (or a continuous transition) is the presence of 
a phase boundary for the former case which separates 
the two phases. The transition for a first order case is 
completed by nucleation of bubbles which expand. The 
moving bubble walls (phase boundaries) lead to physi¬ 
cal phenomena, such as non-trivial scattering of quarks, 
local heating, specific types of fluctuations, etc., which 
are qualitatively different from the case of cross-over or 
a continuous transition. 

It turns out that the presence of moving interfaces is 
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for the case of early universe. 


more generic, and not necessarily restricted to the case of 
first order transitions. Such situations routinely arise in 
the study of so called reaction-diffusion equations 0 ®, 
typically studied in the context of biological systems, e.g. 
population genetics, and chemical systems. Typical solu¬ 
tion of such equations, with appropriate boundary condi¬ 
tions, consists of a traveling front with well defined pro¬ 
file, quite like the profile of the interface in a first order 
transition case. Importance of these traveling fronts in 
the context of high energy physics has been recognized 
relatively recently in several works [5J. In the present 
work we demonstrate such solutions for chiral phase tran¬ 
sition and confinement-deconfinement (C-D) transition 
in QCD even when the underlying transition is a cross¬ 
over or a continuous transition. The only difference be¬ 
tween the field equations in relativistic field theory case 
and the reaction-diffusion case is the absence of second 
order time derivative in the latter case. Thus, correspon¬ 
dence between the two cases is easily established in the 
presence of strong dissipation term leading to a dominant 
first order time derivative term. Further, we will show 
that the required boundary conditions for the existence 
of such a traveling front naturally arise in the context 
of relativistic heavy-ion collision experiments (RHICE). 
For the case of the universe also it may happen in special 
situations as we will discuss below. 

We would like to clarify that the propagating front we 
consider here is like a phase boundary (as in a first order 
transition), and has nothing to do with hydrodynamic 
flow. So, the front will still move, converting one phase 
(say chirally symmetric phase) to the other phase (chiral 
symmetry broken phase) even if the plasma is completely 
static. However, the QGP produced in RHICE undergoes 
hydrodynamic expansion for which one can use either 
Bjorken’s boost invariant scaling model for longitudinal 
expansion, or a 3-d expansion expected to be applicable 
at late stages of plasma evolution. These are incorpo¬ 
rated simply by using appropriate metric for the field 

equations and they lead to a term of the form where 
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the time derivative is with respect to the proper time 
r. It remains to be explored how this new type of phase 
transition dynamics can be incorporated in studies of full 
relativistic hydrodynamical evolution of the plasma. 

There is a wide veracity of reaction-diffusion equations, 
see, e.g. ref. We will discuss specific equations 

which can be identified with the field equations for the 
chiral transition and the C-D transition in QCD in strong 
dissipation limit. Subsequently we will discuss different 
situations in the context of RHICE, with realistic dissipa¬ 
tion, and show that propagating front solutions of these 
equations still persist, making the dynamics of the rele¬ 
vant transitions effectively like a first order phase transi¬ 
tion. 


II. REACTION-DIFFUSION EQUATION FOR 
CHIRAL TRANSITION 

From the form of these reaction-diffusion equations it 
will be clear that such traveling front solutions will exist 
when the underlying potential allows for non-zero order 
parameter in the vacuum state, along with a local max¬ 
imum of the potential Bill- The corresponding values 
of the order parameter provide the required boundary 
conditions for the propagating front solution. First we 
consider the case of spontaneous chiral symmetry break¬ 
ing transition for the two flavor case with the chiral order 
parameter being the 0(4) field </> = (<j,if). We will con¬ 
sider the situation in the context of RHICE, and study 
the transition from chiral symmetry (approximately) re¬ 
stored phase to the chiral symmetry broken phase when 
the partonic system hadronizes during the evolution of 
QGP. For the plasma evolution at this stage, we will 
consider longitudinal expansion as well as spherical ex¬ 
pansion (which may be more appropriate for late stages 
of hadronization). The field equations are 0 : 

<j) — \/ 2 (j) + rj </> = —4A (j) 3 + m{T) 2 (j) + H 
m 2 T 2 

m{T) 2 = -f( 1-^) (1) 

Here, <f> is taken along the a direction. T is the temper¬ 
ature and time derivatives are w.r.t the proper time r. 
We have characterized the dissipation term here in terms 
of 77 which is not a constant for expanding plasma. For 
Bjorken 1-D scaling case rj = 1/t, while for the spheri¬ 
cal expansion 77 = 3/r. We again mention that we are 
only using field equations in the background of expand¬ 
ing plasma where expansion is incorporated by using time 
dependent background metric. It will be interesting to 
study this phase transition dynamics in full relativistic 
hydrodynamical evolution of the plasma. The values of 
different parameters are taken as 6 ], A = 4.5,7710- = 600 
MeV, and T c = 200 MeV. The coefficient of explicit sym¬ 
metry breaking term H = (120MeV) 3 . In the chiral limit 


there is a second order phase transition with the critical 
temperature T c . In view of the explicit symmetry break¬ 
ing term we will take T = 150 MeV which allows for the 
presence of the central maximum in the effective poten¬ 
tial. 

To establish exact correspondence with the reaction- 
diffusion equation, we will first neglect explicit breaking 
of chiral symmetry (i.e. H = 0 in Eq.(l)). Let us also 
consider the extreme dissipative case of large value of 77 
which is time independent, and neglect the cf> term. For 
the resulting equation, we rescale the variables as, x —>■ 
m(T)x , t — y t, and (f> —> The resulting 

equation is, 

0 = V 2 0 - <t? + (t> (2) 

This equation, in one dimension with y 2 </> = d 2 (f>/dx 2 7 
is exactly the same as the reaction-diffusion equation 
known as the Newell-Whitehead equation 0 SI- The 
term d 2 (f>/dx 2 is the diffusion term while the other term 
on the right hand side of Eq.(2) is the so called reac¬ 
tion term (representing reaction of members of biological 
species for the biological systems). 

This equation, in one dimension with y 2 </> = d 2 (f>/dx 2 1 
is exactly the same as the reaction-diffusion equation 
known as the Newell-Whitehead equation 0. The term 
d 2 (f>/dx 2 is the diffusion term while the other term on the 
right hand side of Eq.(2) is the so called reaction term 
(representing reaction of members of biological species 
for the biological systems). We will briefly recall the ana¬ 
lytical traveling front solutions for the Newell-Whitehead 
equation for the present case. Subsequently we will study 
the solutions numerically which will help us in obtaining 
traveling front solutions while retaining the (f> (and with 
77 being time dependent as for the expanding plasma). 

III. PROPAGATING FRONT SOLUTIONS FOR 

CHIRAL TRANSITION 

A. Analytical Solution 

Non-trivial traveling front solutions for the Newell- 
Whitehead equation arise when suitable boundary con¬ 
ditions are imposed, namely (f> = 0 and 1 at x —» ± 00 . 
The analytical solution with these boundary conditions 
has the form, 

cj)(z) = [1 + exp(z/V2)}- 1 (3) 

where z = x — v r. v is the velocity of the front [ 1 ] 
and has the value v = 2>/\f2 for this solution. The 
reaction-diffusion equations typically have several solu¬ 
tions, each with different propagation speeds 0. For 
example, Eq.(2) also has a static solution of the form 
tanh(z). Such a solution can have very important im¬ 
plications for RHICE as well as for cosmology. We will 
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later briefly comment on it. For now we continue with the 
above analytical solution (Eq.(3)). In Fig.la we show the 
propagation of this front. For this we have solved Eq.(2) 
using leapfrog algorithm of second order accuracy. We 
have also added the second order time derivative for nu¬ 
merical solution of Eq.(2) for numerical stability and also 
for comparison with solutions of full Eq.(l). The require¬ 
ment of dissipation dominated dynamics is fulfilled by 
keeping the 77 coefficient of <j> term large, with constant 
77 = 10. This introduces a simple scaling of velocity by 
a factor 1 / 77 . Solid curve in Fig.la shows initial profile 
of the analytical solution in Eq.(3). Plots at subsequent 
times show the propagation of the front. The velocity of 
the front is numerically obtained directly by determining 
the velocity of the front (specifically a particular point 
on the front, say p = 0.5 4> m ax)- We find v = 0.21 for 
77 = 10 in complete agreement with the scaled velocity 

v = - -7=. 

v V2 

For the case of chiral symmetry breaking transition in 
relativistic heavy-ion collisions, the boundary conditions 
required for the solution of Eq.(2) naturally arise due 
to radial profile of energy density of the plasma. The 
center of plasma represents highest temperature T cntr 
which smoothly decreases to values less than the chi¬ 
ral transition temperature T c in the outer regions of the 
plasma. Thus, at any time when T cntr > T c , the chiral 
field will take chirally (approximately) symmetric value 
(which will be zero when H = 0) at center r = 0 and 
will take symmetry broken value at large distances. We 
will take such an initial profile, and evolve it when T cntr 
also reduces to a temperature To below T c . For simplic¬ 
ity, we will assume To to be uniform over the range of 
the profile of <j> with <f> in the center of the plasma hav¬ 
ing a value (f>o (corresponding to the central maximum 
of the potential at T = To). <f> in outer regions of the 
plasma will take vacuum expectation value £ for T = Tq. 
With such boundary conditions, the analytical solution 
in Eq.(3), written in terms of the original (unsealed) field 
and parameters of Eq.(l), takes the form, 

(j>(z)=^[l + exp{ rn ^-(x-VT))]~\ (4) 

where £ = is the vacuum expectation value of p 
(for H = 0) and the velocity becomes v = 3 ^p • Inter¬ 
estingly, the profile of this analytical solution is similar to 
the Wood-Saxon form. From the energy density profile 
expected for colliding heavy nuclei, a Wood-Saxon type 
profile for the field is rather natural. 

It turns out that the form of the traveling front, and 
its evolution, is essentially unaffected even if we take 
non-zero value of H in Eq.(l). Thus, we calculate the 
numerical profile of the front using non-zero value of 
H = (120MeE) 3 . This changes the boundary conditions 
for p(z). For parameter choice in Eq.(l), the vacuum 
expectation value of p is found to be £ = 75.18 MeV 
while the central maximum of the potential is shifted to 


p = <po — -25.93 MeV. The above analytical solution 
in Eq.(4), suitably modified for these changed boundary 
conditions, becomes, 


4>{z) 


(£ ~ M 

Ao 


[1 + exp( 


m(T)(\z\ — Rq) 

V2 


)]- 1 +? (5) 


where the normalization factor Aq = [1 + 

ea :p( ~ rn ^ Ro ]~ 1 . Here we use \z\ in order to have sym¬ 
metric front on both sides of the plasma for the present 1 - 
d case. Rq represents width of central part of the plasma. 


B. Numerical Solution 

We calculate numerical solutions for the full Eq.(l), 
retaining the p term. To correspond to the analytical 
solution, we first consider a large, constant, value of 
77 = (1/0.14)/t 77 _1 in Eq.(l) (even though chiral tran¬ 
sition occurs late at r ~ 4-5 fm), and a uniform fixed 
T = 150 MeV. Note that this does not represent re¬ 
alistic QGP evolution in RHICE. We first study equa¬ 
tions with constant (and uniform) T only to show exact 
correspondence with traditional reaction-diffusion equa¬ 
tions. We will see that the resulting propagating front is 
exactly the same as discussed in literature for reaction- 
diffusion equations. Subsequently we relax this assump¬ 
tion of constant T and study proper time dependence 
of T for expanding QGP. We still retain the assumption 
of uniform temperature for studying front propagation 
as with spatially varying T the effective potential also 
has to vary spatially and correspondence with reaction- 
diffusion equation becomes remote. We consider 1-d case 
as suitable for the Newell-Whitehead equation. This will 
be applicable when size of traveling front is large, so a 
planar approximation can be used. Fig. lb gives the nu¬ 
merical profiles of the front p(z) at different times (with 
initial time taken as 4 fm). The front starts at a dis¬ 
tance of about 10 fm from the center and moves inwards 
converting the central region to chiral symmetry broken 
phase. Note, as To everywhere has a value corresponding 
to the symmetry broken phase, one would have expected 
rapid roll down of the field to p = £ everywhere in time 
of order 1 fm. In complete contrast to this, we see that 
phase conversion here happens slowly, by the motion of 
well defined interface, just like for the case of a first or¬ 
der phase transition. As for Fig.la, we obtain v here also 
directly from the traveling front. We find v ranging from 
0.41 to 0.35 in close agreement with the expected value 
v = (= 0.42). By finding propagating solutions 

for different values of 77 we have verified (for both cases, 
Fig.la, and Fig.lb) that the velocity of the front exactly 
scales as I/77. 

The propagating front solutions we obtain are very ro¬ 
bust and almost independent of the initial profile of the 
front taken. To show this, we show evolution of a pro¬ 
file consisting of linear segments (with correct boundary 



4 



z(fm) 

FIG. 1: (a) Plots of the numerical solution for the traveling 
front of Eq.2 (as discussed in the text) at different times, (b) 
shows the numerical solution for profile in Eq.(5) with non¬ 
zero H. 


conditions) in Fig.2a. We see that this also develops into 
a well defined propagating front as shown in Fig.l. Next, 
we consider realistic values of time dependent r) = 1 /r 
so that second time derivative term becomes important, 
with initial value of r = tq = 4.0 frn, as appropriate for 
the chiral transition. We also take T(t) = Tq^tq/t) 1 ^ in 
accordance with Bjorken’s scaling solution for the longi¬ 
tudinally expanding plasma. Fig.2b shows the traveling 
front solution for this case at different values of proper 
time r. The only minor difference from plots in Fig.l is 
seen at somewhat later stages, with a little rise at the 
boundary of the front. The central value of <f> changes in 
accordance with time dependent T. It is clear that with 
the presence of this front structure most of the results 
for first order transition, such as non-trivial quark scat¬ 
tering, well defined phase separated regions, fluctuations 
etc. will become applicable. 

For Bjorken’s 1-d longitudinally expanding plasma, 
one should consider transverse motion of the front. Ne¬ 
glecting transverse expansion of the plasma, one should 
use Eq.(l) with cylindrical coordinates, or spherical ex¬ 
pansion for the late stages of plasma evolution. Eq. (1) 
for these cases will become 



z(fm) 



z(fm) 

FIG. 2: (a) Initial profile consisting of linear segments also 
rapidly evolves into a well defined propagating front as in (a), 
(b) Solutions of Eq.(l), with realistic values of time dependent 
r / and time dependent T. 


4> ~ 7 T — ~ — 1 ~T + ^ = ~4A</> 3 + m(T) 2 (f> + H (6) 
ar z r dr 

where d = 2 for Bjorken’s 1-d longitudinal expansion, 
and d = 3 for spherical expansion. We have obtained 
front solutions for both these cases here for the chiral 
symmetry case as well as later for the Z(3) case of C- 
D transition. The resulting solutions are very similar to 
those obtained for 1-d case (as in Figs. 1,2). Hence we do 
not show those plots. 

IV. REACTION-DIFFUSION EQUATION FOR 
CONFINEMENT-DECONFINEMENT 
TRANSITION 

We now consider the case of confinement- 
deconfinement transition during early thermalization 
stage. For RHICE, for very early stages, Bjorken scaling 
with longitudinally expanding plasma is a very good 
approximation. The dissipation term is very strong 
during very early stages which helps in making a direct 
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correspondence with reaction-diffusion equations. As the 
traveling wave solutions exist in the symmetry broken 
phase we will consider the case of initial confinement- 
deconfinement transition using the Polyakov loop order 
parameter with center symmetry Z(3) spontaneously 
broken in the high temperature QGP phase. This early 
stage actually represents a non-equilibrium stage, with 
system thermalizing to a maximum temperature To i n 
a time scale tq. The tliermalization time scale tq can 
be as short as about 0.14 fm (at LHC). Elliptic flow 
measurements indicate an upper bound of about 1 fm 
on To- For the dynamics of the order parameter l(x) 
which is the expectation value of the Polyakov loop, we 
use the following effective Lagrangian density [7j . 

L = ^|Zy| 2 T 2 - V(l) (7) 

Where the effective potential V(l) for the Polyakov 
loop, in case of pure gauge theory is given as 

v(i ) = (^|/| 2 - |a 3 + (o 3 ) + \m 2 ) 2 )h t 4 (8) 

At low temperature where l = 0, the potential has only 
one minimum. As temperature becomes higher than T c 
the Polyakov loop develops a non vanishing vacuum ex¬ 
pectation value Z 0 , and the l 3 + l* 3 term above leads to 
Z( 3) generate vacua. Now in the deconfined phase, for 
a small range of temperature above T c , the l = 0 ex¬ 
tremum becomes the local minimum (false vacuum) and 
a potential barrier exist between the local minimum and 
global minimum (true vacuum) of the potential. As we 
are interested in showing the existence of traveling front 
solutions in the absence of any first order transition, we 
will consider value of temperature T to be sufficiently 
large so that there is no such barrier present. For the pa¬ 
rameter values we use, this requires T > 280 MeV, and 
we take T = 500 MeV. The values of various coefficients 
in Eq.( 8 ) are the same as used in our previous works [5112] 
(including discussions about explicit symmetry breaking 
strength &i) and we do not repeat that discussion here. 
(With those values of parameters, the transition temper¬ 
ature is taken to be T c = 182 MeV.) 

For simplicity we neglect the effect of dynamical quarks 
which lead to explicit breaking of Z(3) symmetry, and 
hence a linear term in l in V(l) above [JJ]. This can be 
taken care of in a similar manner as the explicit symmetry 
breaking term H for the chiral symmetry case in Eq.(l). 
Similarly, as our interest is not in the Z(3) structure of the 
vacuum, we will take l to be real. We again first neglect 
the second order time derivative (for large dissipation 

case). The variables are scaled as, x — > gTx , and 

7 2 m2 

t —> *9 N t. With that, the field equation for (real) l(x) 
can be written as follows (for the sake of uniformity, we 
will denote l(x) as <j>(x) in the following), 


</> = V 2 ^ + 4>( b 2 + b 3 (j) - </> 2 ) (9) 


The final equation in this case is again a reaction- 
diffusion equation known as the Fitzhugh-Nagumo equa¬ 
tion which is used in population genetics El H|. Thus 
we again expect well defined traveling wave solutions for 
appropriate boundary conditions. 


V. PROPAGATING FRONT SOLUTIONS FOR 
CONFINEMENT-DECONFINEMENT 
TRANSITION 

Required boundary conditions for the propagating 
front solution for Eq.(9) again naturally arise in RHICE, 
during early stages. As the system thermalizes, one ex¬ 
pects first the center of the plasma to reach a tempera¬ 
ture T > T c where T c is the C-D transition temperature. 
The temperature in (somewhat) outer regions remains 
below T c initially. This leads to a profile of l(x) where 
l = 0 in the other regions while l = 1 at the center 
of the plasma. Subsequently even these regions, some¬ 
what away from the center, also achieve T > T c . With 
these boundary conditions, we solve equations for l(x) 
with a uniform temperature T with initial value = 500 
MeV. Eq.(9) above was derived in the large dissipation 
limit to identify it with the Fitzhugh-Nagumo equation 
which guarantees the existence of a traveling wave so¬ 
lution for l{x). With that assurance, we will now solve 
the full field equations for l(x) (i.e. for (f>(x)) including 
the second time derivative term. The dissipation term 
is naturally large initially in this case due to 1/r factor 
with tq being very small. Again, to show direct corre¬ 
spondence with reaction-diffusion equation, we will take 
a very large, fixed, value of 77 = 1 /tq with Tq = 0.01 fm, 
and will keep temperature to have fixed value T = Tq = 
500 MeV. 

Plots in Fig.3a show the well defined traveling wave 
solution at different values of r starting from initial time 
taken as To = 0.14 fm. The initial profile is taken to have 
similar form as in Fig.l, suitably modified for the bound¬ 
ary conditions appropriate for the present case. We next 
consider realistic value of time dependent ?/ = 1/r with 
initial value of r = To = 0.14 fm, and take T(r) = 
ToIto/t) 1 / 3 as appropriate for the Bjorken 1-d scaling 
solution. Resulting evolution of <j>(z) is shown in Fig.3b. 
Though (f> shows some oscillations, still it shows a reason¬ 
ably well defined propagating front. It is possible that rj 
may not decrease as fast as 1 /r due to presence of other 
sources of dissipation. In that case resulting solution will 
be closer to that in Fig.3a. 
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FIG. 3: (a) Traveling wave solution for the Polyakov loop 
order parameter at different times for very large dissipation 
case with constant ?? = 1/(0.Olfm). (b) Solution for realistic 
rj = 1/t and with time dependent T. 


VI. RELATION BETWEEN CHIRAL 
TRANSITION AND DECONFINEMENT 
TRANSITION 


We now address an important issue regarding the dif¬ 
ference in the behavior of order parameter evolutions for 
the chiral symmetry case and for the Polyakov loop for 
the C-D transition case. Note that we discussed trav¬ 
eling wave solution for the chiral symmetry case during 
late stages of the evolution, even though presence of large 
dissipation case was harder to justify for such late stages 
(1/r factor being relatively smaller). Similarly, for the 
early thermalization stage, we only discussed the case of 
C-D transition with the Polyakov loop, and did not dis¬ 
cuss the chiral symmetry case. The reason is that the 
traveling front solutions via the reaction-diffusion equa¬ 
tion approach arises only when the effective potential has 
specific shape, for example, that corresponding to spon¬ 
taneous symmetry breaking. Thus, spontaneous break¬ 
ing of Z(3) symmetry during early thermalization leads 
to well defined traveling front for the Polyakov loop. But 
during this stage, chiral symmetry is restored. There 


is absolutely no possibility of finding any traveling wave 
solution for the symmetry restored effective potential for 
the chiral field </> (as one can simply check analytically 
as well as numerically). Similarly, while a traveling front 
exists for the chiral field during chiral symmetry break¬ 
ing transition during late stages of the plasma, there is 
no such solution for the Polyakov loop order parame¬ 
ter at that stage as the Z(3) symmetry actually gets re¬ 
stored. This raises serious concerns with the conventional 
idea that the two transitions, namely the chiral transi¬ 
tion and the C-D transition, are somehow related (or, 
are the same). In our calculations we clearly find regions 
separated by the traveling front of one order parameter 
whereas no such phase separation is expected from the 
other order parameter. Thus, we will conclude that our 
results support the claims of several groups that the chi¬ 
ral transition and the C-D transition are indeed separate, 
and independent, transitions. In fact, our results provide 
clearly phase separated regions (by the traveling front) 
which are, e.g. chirally symmetric but in confined phase, 
or in chiral symmetry broken phase but in the deconfined 
phase. 


VII. CONCLUSIONS 

We conclude by emphasizing that the techniques of 
reaction-diffusion equation have been used here to show 
existence of well defined traveling front solutions, which 
are very similar to phase boundaries for a first order 
transition case, even though the relevant QCD transi¬ 
tions here are of second order, or a cross-over. This 
allows the very exciting possibility of using earlier re¬ 
sults valid for a first order transition case, such as for¬ 
mation of strangelets, baryon concentration, fluctuations 
etc. for RHICE. Our results show that the transition 
proceeds by slow moving front, and may take several fin 
time to complete, leading to long lasting mixed phase 
stage. This will affect calculations of various signals of 
QGP for RHICE, e.g. production of thermal photons 
and di-leptons, J/ip suppression, and especially elliptic 
flow which develops mostly during early stages. As we 
mentioned above, the reaction-diffusion equations have 
other solutions also with different propagation speeds. 
For example, Eq.(2) also has a static solution of the form 
tanh(z). If the initial density profile in RHICE leads to 
such a profile of 4>{z) the transition will become stag¬ 
nant. This will have important implications for RHICE. 
It is clearly of great importance to see if these results 
can be applied to the case of early Universe (for exam¬ 
ple, solutions with different speeds, especially the one 
with zero speed). The required initial boundary condi¬ 
tion of (f> = 0,1 for x —±oo looks difficult to justify for 
the Universe (due to absence of a temperature profile as 
for RHICE). However, one should remember that there 
are always density (and hence temperature) fluctuations 
present in the Universe (most likely of inflationary ori¬ 
gin). Though these are very tiny (one part in 10 5 ), but 
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if we consider these fluctuations when the temperature 
of the Universe is very close to the transition tempera¬ 
ture T c then there can easily be regions of space where 
the symmetry is restored, while other neighboring regions 
will have symmetry broken phase. 

This will lead to the required boundary conditions for 
the traveling wave fronts as discussed here. Small magni¬ 
tude of temperature fluctuations will imply small differ¬ 
ence in the magnitude of $ at the two boundary points, 
subsequently leading to small effects (like scattering of 
quarks). However, typical wavelength of these fluctua¬ 
tions will be comparable to the Hubble size, naturally 
leading to wavefront propagation over such large scales. 
Thus, even with small quark scattering etc. one may be 
able to get large concentration of baryons via Witten’s 
mechanism of nugget formation. 

One caveat in this scenario is that chiral symmetry 
(as well as Z(3) symmetry) are also explicitly broken. 
This makes the above mentioned scenario difficult to im¬ 
plement for the quark-hadron transition when tempera¬ 
ture fluctuations have very small magnitude. There may 


be other possibilities for quark-hadron transition which 
we intend to explore in a future work. However, such 
a mechanism will certainly apply for electroweak phase 
transition where there is no explicit symmetry breaking 
involved. There also one will get traveling front solu¬ 
tions arising from inflationary fluctuations irrespective 
of the order, or the strength of phase transition. It will 
be interesting to investigate how such front solutions can 
affect the physics of post electroweak transition physics, 
in particular sphaleron mediated baryogenesis etc. 
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